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^ ' Abstract 

a: 

q^' Quarkonium production at a given large pr is dominated by parton 

D ! fragmentation: a parton which is produced with transverse momen- 

tum pt/z fragments into a quarkonium state which carries a fraction 
z of the parton momentum. Since parton production cross sections 
^ , fall steeply with pt, high z fragmentation is favored. However, quan- 

tum number constraints may require the emission of gluons in the 
fragmentation process, and this softens the z dependence of the frag- 
mentation function. We discuss the possibility that higher-order pro- 
cesses may enhance the large z part of fragmentation functions and 
thus contribute significantly to the quarkonium cross section. An ex- 
plicit calculation of light quark fragmentation into r\ c shows that the 
higher-order process q — » qr\ c in fact dominates the lowest-order pro- 
cess q ->■ qgr] c . 



1. Introduction 

Heavy quark production is a hard QCD process, and the perturbative ex- 
pansion of production amplitudes in a s is expected to apply. However, the 
perturbation series is not always dominated by its lowest-order term. For 
example, when the transverse momentum of a collinear heavy quark pair 
is much larger than its invariant mass M, light parton fragmentation dom- 
inates the lowest-order production processes by a power oip\jM 2 HjJ. The 
data on charmonium production indeed shows a l/p? behavior of the cross 
section @, ||, compatible with the prediction for a fragmentation process. 

QCD calculations based on the color singlet model [f| nevertheless dis- 
agree with data on the relative production rates of S and P wave quarkonia 
and on the absolute normalization of the cross sections at both low || |6| 
and high ^, ^] values of px- Order-of- magnitude discrepancies have been 
observed for both charmonium and bottomonium. 

In QCD, the fragmentation of a virtual gluon into a 3 Si quarkonium state 
(J/ip, T) requires the emission of at least two extra gluons, g* — > 3 Si + gg, 
whereas fragmentation into a 3 Pj state can proceed with the emission of a 
single gluon, g* — > 3 Pj + g. The emission of the extra gluon suppresses the 
calculated 3 Si cross section considerably compared to the 3 Pj cross section, 
due to the extra power of a s and also because the emitted gluons carry away 
part of the transverse momentum of the fragmenting gluon. The experimen- 
tal 3 Si/ 3 P2 cross section ratio is much larger than the calculated one. 

The discrepancy could be related to the bound state dynamics of the 
quarkonia. This is the solution proposed by the color octet model [§, [|. The 
gluon is assumed to first fragment into a QQ pair in a color octet state. Later, 
after a formation time characteristic of the bound state, the pair couples to 
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the physical quarkonium through the absorption or emission of one or two soft 
gluons. The momenta of the emitted gluons are typical of the quarkonium 
bound state dynamics, and thus they carry away only a minor part of the 
transverse momentum. 

Although the probabilities of these nonperturbative transitions are essen- 
tially free parameters, the octet model makes specific predictions about the 
polarization of the produced quarkonium flIU| . The 3 Si quarkonia produced 
by the fragmentation of a nearly real gluon are expected to be transversely 
polarized. The data to test this prediction is not yet available. Note that the- 
oretical calculations of quarkonium polarization can be compared with data 
in fixed-target experiments. In this case, neither the color singlet model nor 
the color octet model can explain the experimentally observed unpolarized 
production |], [HJ . 

Here we wish to draw attention to the possibility that higher-order per- 
turbative contributions to fragmentation mechanisms could be enhanced by 
the so-called trigger bias effect. When the transverse momentum pt of the 
quarkonium is fixed, processes which allow the fragmenting parton to be pro- 
duced with the lowest possible transverse momentum Pt/z are favored. The 
large px cross section is a convolution of the production cross section <Ji of a 
parton i and its fragmentation function D^q to the quarkonium state O, 



dcr (s,p T ) _ r l ^ da, 
dp T ^ Jo dp iT 



I dz-^-(s,p T /z,fx)Di^ (z,n), (1) 

: Jo dVirn 



where [i is the factorization scale. The parton production cross section 
d<Ji/dpi T falls approximately as (pt/z)~ a , which implies the enhancement 
of high z fragmentation by a factor z 4 . As a rough measure of the impor- 
tance of a fragmentation function Di^o(z, /i) we can therefore use its fifth 
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moment D^ (fi), denned by 

f dzz^D^oiz^). (2) 
Jo 

We have studied the importance of higher-order fragmentation contribu- 
tions in the case of light quark fragmention into 1 So quarkonium (i] c ) within 
the color singlet model. The relevant Feynman diagrams are shown in Fig. 
[I| In the nonrelativistic limit the fragmentation function is a product of the 
fragmentation probability into a collinear, on-shell cc pair in a 1 Sq state and 
the square of the r\ c wave function at r = 0, 




g 

Figure 1: (a) A lowest-order diagram contributing to the process q — > 1] C +X. 
(b) A higher-order diagram with no gluon emission. In each case there is 
another diagram with the c-quark-gluon vertices interchanged. 

At lowest order, q — > r] c + X fragmentation is due to the process q — > qr] c g 
of Fig. |l|a. The emission of the gluon suggests that this process may have a 
softer fragmentation function than the higher-order process q — > qr] c shown 
in Fig. |I|b. Due to the trigger bias effect, the higher-order process could be 
enhanced. 
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2. Results 



Our calculation of the q — > rj c + X fragmentation processes shown in Fig. |l| 
is described in the Appendix. The contribution of the lowest-order process 
(Fig. |l|a) to the fragmentation function is of the form 



The coefficient functions are 

(6 (2 + z) 



f(z) 



a 3 s C F \R s (0) t 
48n 2 m 3 N r 



4m 2 c 



7l~ 



4m 2 



(4) 
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+ --18 + 12z + 4z 2 + 6 + --12z ln(l - z 



C 2 (z) 



6 



3 z ha.(z) 



(5) 



9{z) 



it 2 (2 + z) ln(z) + -18 + - + 18 z + 4 z 2 ) ln(z 



a 3 s C F \R s (0)t 
48ir 2 m 3 c N c 

+ U+--12z\ ln(l-z) ln(z) - 3 ^ln^) 2 + 34 + vr 2 (- - 2 

71 53 ^ 13 z 2 /-5 n . 2 \ . H 
-6-z- — + — + {- + 9z - AZ ) H1 ~ Z) 

+ (18-12z) C 2 (z) + 6 (2 + z) [£ 3 (z)-\n(z)C 2 (z)-((3)}\, (6) 



where 



C 2 (z) 



ln(l - t) 



dt 



is the dilogarithmic function, 



£a(z) 



o t 



(7) 



(8) 



and C(3) « 1.202. The logarithmic term /(z) ln[/i 2 /(4m 2 )] arises from the 
two-step process where q —>■ qg splitting is followed by g —>■ rj c g fragmentation; 
the function f(z) can be written as 

dy 



/(*) 



* y 



■P^(z/y)D^Jy), 



(9) 



where P q _ g {z/y) is the standard q — > qg splitting function [13] and D (y) is 
the g ?7 c f7 fragmentation function at lowest order M. A similar result has 
been obtained in the case of J/ip production by light quark fragmentation 
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A lower limit of the loop contribution (see Fig. is obtained by consid- 
ering only the imaginary part of the loop amplitude. There is no logarithmic 
term in this case: 



where 



h(z) 



D^Jz,fi)>h(z) + 
ai\R s (0)\ 2 C 2 F 



96nm 3 N r 



+ z + 



2z 



ln(z) + 



14(1 -z) 

z(7z 2 



4m^ 
6 



(10) 



Cod 



18z + 12) 



W(z) 



11) 



l-z~" y '' ' (l-z) 2 , 
The functions f(z), —g(z) and h(z) are plotted in Fig. using a s = 0.26, 
|-Rs(0)| 2 = (0.8 GeV) 3 , and m c = 1.5 GeV. The loop contribution dominates 
over the lower-order Born contribution for z <; 0.3, even though the real part 
of the loop was neglected. More quantitatively, the fifth moments of the Born 
and loop contributions have the numerical values 



J 

DfX^)> f Q dz*h{z) 



2.4 In 



-5.1 



x 10 



-7 



1.1 x 10" 



(12) 
(13) 



Depending on the fragmentation scale /i, the contribution from the loop 
diagram is thus up to an order of magnitude larger than the lowest-order 
Born contribution. Neglecting the higher-order process would lead to a major 
underestimate of the fragmentation cross section. 

It is possible to further simplify the calculation of the fragmentation func- 
tions by taking advantage of the fact that only the large z region is important, 
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Figure 2: The functions f(z), —g(z) and h(z) as defined in the text. 
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due to the trigger bias effect. We have verified that using only the leading 
part of an expansion of D(z,fi) around z = 1 changes the fifth moments of 
the loop and Born contributions by less than 10%. 

3. Discussion 

The trigger bias effect in large px quarkonium production favors fragmenta- 
tion processes where the quarkonium takes a large fraction z of the momen- 
tum of the fragmenting parton. When estimating the relative importance 
of different fragmentation processes, the shape of their z dependence must 
therefore be considered. 

In particular, some higher-order perturbative contributions may be en- 
hanced relative to the lowest-order contributions due to the trigger bias effect. 
In this paper, we analyzed the process q — > r\ c + X, where such an enhance- 
ment can be expected because gluon emission is not required in higher-order 
processes. We found that there is a loop contribution which indeed dominates 
the Born contribution by a large factor. 

It is likely that an analogous result is obtained in the case of q — > J/ip 
fragmentation. Some relevant Born and loop diagrams are shown in Fig. 
|3]. At higher orders, all the gluons coupling to the heavy quark line can be 
attached to the light quark line instead of being emitted, which suggests a 
hard z dependence of the fragmentation function. 

These higher-order contributions are part of the standard perturbation 
series and thus do not bring in any new parameters. Their relative impor- 
tance should depend only weakly on the quark mass (through the decrease of 
Q^O^q) with tuq). This is in qualitative agreement with total cross section 
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> JA|f 



a 



Figure 3: Light quark fragmentation into a J/i/j. 

data [fj, [| ||, which shows a disagreement with Born term calculations 
(within the colour singlet model) of similar magnitude for bottomonium and 
for charmonium. 

The calculation presented here is not, however, immediately applicable 
to the present data on quarkonium production. The primary production 
mechanism for quarkonia at large pt in hadron collisions is expected to be 
gluon fragmentation. Even at higher orders, a minimum of two extra gluons 
need to accompany a produced J/ip, due to charge conjugation invariance (c/. 
Fig. D. In this case, loop diagrams like the one in Fig. [|b simply represent 
radiative corrections to the lowest-order process. Whether they enhance the 
kinematic region where the emitted gluons carry little momentum (the large 
z region) can only be determined by an explicit calculation. 

On the other hand, processes such as the one in Fig. |3|b could be signifi- 
cant in collisions where light quarks are more copiously produced relative to 
gluons, such as at HERA. There, however, also charm quark fragmentation 
becomes important as a charmonium production mechanism at large pr ]15 . 
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In summary, we have pointed out that the trigger bias enhancement of 
large z fragmentation is crucial in quarkonium production at large p?- As 
a specific example, we considered the q — > rj c fragmentation process and 
calculated a higher-order perturbative correction whose contribution to the 
cross section exceeds the lowest-order fragmentation contribution by a large 
factor. 

Acknowledgement. We are grateful for discussions with Stan Brodsky. 




Figure 4: Gluon fragmentation into a J/ip. (a) A lowest-order diagram, (b) 
A higher order diagram. 



Appendix 

We describe here our calculation of the q —*■ rj c fragmentation functions. As 
shown in Figs. |5| and |6]a, we denote by p the momentum of the quarkonium 
state; Q denotes the momentum of the fragmenting light quark, and s = Q 2 
is its virtuality. 
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We work in the center of mass system of the light quark production 
process, and choose the third axis along the direction of the fragmenting 
light quark. The light-cone components of a four-vector v are defined as 
v = v° ± v 3 and v± = (v l ,v 2 ). The variable z is defined as z = p + /Q + , 
which in the limit of large Q + is the fraction of the light quark momentum 
taken by the r\ c . 

We use an axial gauge with the polarization tensor 

h^n u -A- n^h v 

V» v (k,n)=g» v - + , (14) 

n ■ k 

where the gauge vector n satisfies n 2 = 0, and n ■ v = v + /Q + . 

Let us write the matrix element for light quark production as u a (Q)A4 a , 
where a is a Dirac index. The square of the amplitude for the full r] c pro- 
duction process can then be written as A4*^Tp a A4 a . In the limit of large 

T fBa = T$ Pa + ... (15) 

where T is a scalar function. The full cross section then becomes a convolu- 
tion of the light quark production rate 

1 

V CM 

and a q — *■ r\ c fragmentation function which is given by a phase space integral 
of T, as shown below. 



da = TTE^- dLi P s M *p$ p a M c 



(16) 

Q 2 =0 



The leading-order fragmentation function 

At lowest order, the q — > rj c fragmentation function gets contributions only 
from the Feynman diagram of Fig. ^ and another diagram where the two c- 
quark-gluon vertices have been interchanged. The amplitudes corresponding 
to the two diagrams are equal. 
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As shown in Fig. |5], we denote the momentum of the virtual gluon by k 
and define w = k 2 and y = n ■ k. 




Figure 5: Momentum definitions in light quark fragmentation into f] c . 

Let us first consider the square of the q* qrj c amplitude as it appears 
in the cross section of the full process. We find 



T 



1287r 2 a 3 s C F \R s (0)\ 2 



m c N c 



-JO 



x 



V aa '(k,n) ppk^e 



k 2 k 2 — Am 2 
T Q + T p] / + T k H + T n ri 

(Tq + zT p + yT k ) l -Q + T + ... 
{T Q + zT p + yT k ) Q + . . . 



s 



k 2 — Am 2 



k 2 



(17) 



The tensor p Q £^e Q/3/ ^ is due to the cc spin projection [|l6j. The dots in the 
last two expressions stand for terms of relative order 1/Q + . We made use of 
the fact that the coefficients Tg,T p ,T k , T n depend only on scalar products of 
the four- momenta and are therefore independent of Q + . Explicitly, 



T Q + zT p + yT k 
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128n 2 aiC F \Rc 



C- 1 ' c 

1 



X 



2s 2 w 2 (w — Am 2 ) 2 



-32 m\ s — 16 m\ w — A pi s w — 2 s w z 



+ 



2 s (lQmf + 2plw + w 2 

y 



+ s (l6 m\ + w 2 ^j y 



2 (A m 2 c +pl) 2 w 



2 (Ami + pi) w (Ami + w) 2 sw (Am 2 + w) z 2 2 
I — s w z 



y 



2 s 2 w z 2 ^2sw (2 s + w) z 2 



(18) 



v y 

The phase space measure for the full process can be written as the product 

of three factors: the phase space measure for light quark production, the 

phase space measure for the decay of the virtual light quark, and ds/(27i). 

We write the two latter factors as 
ds d 4 k d A p 



(2tt) (2tt) 4 (2tt) 



4 2ttS + ((Q - k) 2 )2v5 + (p 2 - Am 2 c )2n5 + ((k - pf)9(p - s) 



dz 



1 



ds — 



dy f s v 



dw 



dt 



2567T 4 Hmllz J z y Jim? c y/z J-Jp TTy/p — t' 2 J (27r) 

where <fi is the azimuthal angle of p±, and 

Az 2 



P 
t 



y 



■(1 - y){y - z)(sy - w) (wz-Am 2 c y) 



w z 



(2z — y — y z) H — (s z — Am c ) + Am c + s z 



(19) 



(20) 



(21) 



y y 

The integral over z gives the convolution in the production cross section, and 
the leading-order light quark fragmentation function is 

2567T 4 JAml/z Jz y JAmly/z J-y/p 7Ty p — t 

m^(^)+9(z) + 0(^]+0(at). (22) 



Am 2 



Analytical expressions for the functions f(z) and g(z) are given in eqs. 
and (El), respectively. 
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The loop contribution 

The trigger bias enhanced NLO contribution to the q —>■ r] c fragmentation 
function comes from the Feynman diagram in Fig. |6]a, and another diagram 
where the two c-quark — gluon vertices have been interchanged. The ampli- 
tudes corresponding to these two diagrams are equal. 

The four-momenta are defined in Fig. |5]a. The loop momentum is denoted 
by k, and y = n ■ k. 




p/2 

— w 



QS k -p/2 




a 



Figure 6: (a) Momentum definitions for the loop diagram contribution to r] c 
production, (b) The cuts which give the imaginary part of the loop diagram. 
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We first consider the structure of the box loop integral 

d A k 7^$ -^7° V afi {k,n) V Pv {p-k,n) p a k p e a ^ 



x 



x 



• (23) 



(2tt) 4 (Q-k) 2 k 2 (p-k) 2 (k 2 -p-k)' 

The four factors in the integrand of eq. (|23|) are easily identified with the 
four sides of the box loop of Fig. |^. Making a Dirac decomposition of the 
integrand we find 

d 4 k 



(24) 
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where 



den x a M 



-V (k-p) - z(k-Q) + 2m 2 y + — + k 2 z k a n p p^ 
- z(k-e) k^ + y(k-e)p^, 
(1-y) Wa\y + (z - 2y)(k-p)] k^ 
-(l-y)\y(k-p)- k\2y-z)\p» 



.a/3-yfj, 



(25) 



den x If 




and the denominators from the propagators and the gluon polarization ten- 
sors are included in 



Any four- vector can be written as a linear combination X M = X n n^ + 
XqQ^ + X p p^ + X e e M of n, Q, p and e M = n a Q pp^e"^ 11 . It is easily seen that 
a p (k), aq(k), a n (k), and b e (k) are all antisymmetric when k is mirrored in 
the hyper plane spanned by n,p and Q, i.e. when k e — > —k e . Therefore they 
do not contribute to the integral, and 



den = y{z - y)(Q - k) 2 k 2 {p - kf k 2 - (p-fc) 



(27) 




(28) 



B» = B n n» + B Q Q» + B p p» 

= / -Si M*K + b Q (k)Q» + b p (k)p»] . 



(29) 



Analogously with eq. (|17D , we now find 
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K Q + K p f + K n f + iKj^ff 
(K Q + zK p )0 + ... 

8iv c m£ s z 



(30) 



where 



C = (sz-Am 2 c )A e -(2B n + sB Q + sB p ) 
d 4 k 

——j (sz - Am 2 c )a e - (2b n + sb Q + sb p ) . (31) 
(2tt) 4 l • 

The phase space measure for the decay of the virtual light quark times 
ds/(2n) is in this case 

1 /■! rU 2 /-27T 



, , dz I ds 

327T d JO Jiml/z JO 



(32) 



As C is independent of the azimuthal angle of the decay, we obtain the 
following 'box' contribution to the light quark fragmentation function: 

:i-z) 



Dibox)(z aj\Rsm Z Cl f 

1287ri V cm 5 J 4m 2 /z - s 2 



ds- 



■\c\\ 



(33) 



We have only calculated the imaginary part of the box amplitude, which 
is due to the sum of the two Cutkosky cuts I and 77 of Fig. This gives a 
lower limit of the full loop contribution. The imaginary part is obtained by 
replacing, respectively, 

d 4 k 



d 4 k 



^ 1 

~* 2 

ii 1 

(2tt) 4 ~* 2 



(2vr) 4 
d 4 k 



U-^ k2{p ~ kf 2n6+{k2) 2n5+i(p - k)2) 



(2vr) 4 
d 4 k 

W) 



k 2 [Q - k) 2 27r6 + (k 2 ) 2u5 + {{Q - k) 2 ) (34) 



m eq. 



ID. We find 
Im(C) = C T + C n 



C 1 = —In 



Am c s 
s — Am 2 
fs — Am 2 )z 



2s - Am 2 , 

1 T^r ln 

s — Ami 



Am 2 



s z 



Am 2 



(35) 
(36) 
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Performing the s integration we get the result given in eqs. ( 
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